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Abstract. PHENIX measurement of the fractional momentum loss (iIpt/pt) of high px 
identified hadrons are presented. The ^px/px of high px rr® which are computed from 39 GeV 
Au-|-Au over to 2.76 TeV Pb-|-Pb are found to vary by a factor of six. We plotted the 5px/px 
against several global variables, Apart, Aqp and dNch/dr], and found global features. It was 
found that 200 GeV Au-I-Au points are merging into the central 2.76 TeV Pb-fPb points when 
plotting Spt/pt against dN^h/dp. 


1. Introduction 

The interaction of hard scattered partons with the medium created by heavy ion collisions (i.e., 
quark-gluon plasma, QGP) has been of interest since the beginning of the RHIC running [T]. 
A large suppression of the yields of high transverse momentum (px) hadrons which are the 
fragments of such partons was observed, suggesting that the matter is sufficiently dense to cause 
parton-energy loss prior to hadronization [2]. The PHENIX experiment [3] has been exploring 
the highest px region with single vr® mesons which are leading hadrons of jets. We show a 
calculation on the energy loss of partons published almost 20 years ago in Fig. [D^a). Although 
the measurement of the momentum shift is the ultimate goal, the paper suggested looking at 
the ratio of the high px single hadrons in Au-|-Au and p+p collisions as an alternate way. 
Since then, most of the experiments including PHENIX have looked at the nuclear modification 
factors, Raa (= {dNAA/dydpT)/{{TAA)dcrpp/dydpT)), and quantihed the energy loss effect via 
its suppression. We here present the momentum shift of high px hadrons instead of Raa- 

2. Fractional momentum loss dp^/pT 

With a larger statistics of both p -|- p and Au-|-Au data recently collected, it became possible 
to measure the momentum shift directly. Fig. [U^b) depicts the method to compute such shift. 
We have statistically extracted the fractional momentum loss {6pt/pt, dpT = px — Pt', where 
Px is the transverse momentum of the p -|- p data, and px’ is that of the Au-|-Au data) of the 
partons using the hadron px spectra measured in p -|-p and Au-|-Au collisions [3]. Since the 
number of data points is finite, a fit to the scaled p -|-p is needed to evaluate dpx/pT at a given 
Au-|-Au invariant yield. The uncertainty of the 5px/PT is calculated by inversely converting the 
quadratic sum of the uncertainties on the yields of Au-|-Au and p -|- p points, using the p -|- p fit 
function. Statistical and systematic uncertainties are individually calculated in the same way. 


1 




Figure 1. (a, left) A calculation demonstrating that the suppression of the vr^ yield is resulted 
in from the shift of the momentum spectra in A+A collisions, (b, right) Method of calculating 
average (5 pt/pt- We scaled the p + p yield by Taa corresponding to centrality selection of 
Au+Au data, shifted the p + p points closest to Au+Au in yield, and calculated the momentum 
difference of p + p and Au+Au points. 


The uncertainties on Taa and p+p luminosity are not plotted but mentioned in plots. Fig. [JJ^a) 
show the Raa for the vr^’s in Au+Au collisions at =200 GeV from the RHIC Year-7 run. 

Using the Sp^/pT calculation method, we obtained the 6pt/pt for the same dataset as shown 
in Fig. [2Db). Similarly, the Raa for the tt^’s in 0-10% Au+Au collisions at =39, 62 and 
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Figure 2. (a, left) Raa of vr^’s for 200 GeV Au+Au collisions obtained from RHIC Year-7 run. 
(b, right) 5pt/pt for the same dataset. 




















200 GeV from the RHIC Year-7 and Year-10 runs shown in Fig. [3]^ a) are replotted in the form 
of 5pt/pt as shown in Fig. [3)^b). The Raa^ look similar even the cms energies are changed 
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Figure 3. (a, left) Raa of vr'^’s for 39, 62 in 0-10 % Au-|-Au collisions obtained from RHIC Year- 
10 run and in 0-10% 200 GeV Au-|-Au collisions from RHIC Year-10 run. (b, right) 6pt/pa for 
the same dataset. 


by a factor of 5 as seen in Fig. [3Da). However, when we look at SpaIpt for the corresponding 
dataset, we found that the Sp^/pA changes by a factor of three from 39 to 200 GeV as shown in 

Fig. EDb) [5]. 

The Raa's also look similar between RHIC and LHC (Fig. HKa)). Similarly, the Spa/pa is 
found to change by a factor of ~ 1.5 from 200 to 2.76 TeV (Fig. HKb)). To summarize, even the 
Raa^s are similar, the (5 pt/pt’s show a factor of six variation from 39 GeV to 2.76 TeV. This 
fact has not been found by looking at Raa- 


3. Scaling property of 5pa/pa 

In order to study the systematics of 6pa/pa, we plot the Spa/Pa against several global variables 
such as Vpart) A^qp (number of quark participants) and dNch/dp. We first plotted the Spa/pa 
against Apart as shown in Fig. [5l All the plots shown in this section are at pa{p + p) = 7 GeV/c 
in order to reach the hard scattering regime. In the Apart scaling we see that the Cu-|-Cu and 
Au-|-Au are nicely lined up, implying that within the same cms energy, the Spa/Pa scales with 
Apart- This is consistent with the fact that Raa is similar at same Apart between Cu-|-Cu and 
Au-|-Au collisions [7]. The Pb-|-Pb points are consistently off the trend of 200 GeV points, but the 
slopes of both systems look similar. Fig. [6] shows Spa/pa against Aqp. The detail description of 
how the number of quark participants are obtained can be found in the literature [8]. We 
employed a Monte-Carlo-Glauber (MC-Glauber) model to calculate the numbers. We first 
determine the quark-quark inelastic cross section (cr“®*) for each collision energy such that 
the inelastic nucleon-nucleon cross section {up/fj) is reproduced. Then the model is modified to 
handle the quark-quark rather than nucleon-nucleon collisions. The nuclei are placed according 
to a Woods-Saxon distribution and then three quarks are distributed around the center of each 










Figure 4. (a, left) Raa of vr^’s for 200 GeV Au+Au collisions obtained from RHIC Year-7 
rnn and charged hadrons for 2.76 TeV Pb+Pb collisions obtained by the ALICE experiment at 
LHC [6]. (b, right) 5pa/pt for the same dataset. 



• Pb+Pb 2.76TeV, Ch. hadrons: PLB 720, 52 (2013), PRL 106, 032301 (2011) 


• Au+Au 200GeV,n°: PRC 87, 034911(2013) 

• Cu+Cu 200GeV,n°: PRL 101, 162301(2008) 
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Figure 5. (Ipt/pt as a function of Apart for vr^’s in 200 GeV Au+Au and Cu+Cu collisions 
measured by PHENIX and charged hadrons in 2.76 TeV Pb+Pb collisions measured by ALICE. 


nucleon following the distribution of: 

pproton^^^ = X 6-“’’ 

where a = y/vlirm = 4.27fm~^ and r„i = 0.81 fm is the rms charge radius of the proton. A 
pair of quarks, one from each nucleus, interact with each other if their distance d in the plane 
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• Pb+Pb 2.76TeV, Ch. hadrons: PLB 720, 52 (2013), N calculated by PHENIX 

• Au+Au 200GeV,7i“: PRC 87, 034911(2013) 

• Cu+Cu 200GeV,7i“: PRL 101,162301(2008) 
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Figure 6. 5pt/pt as a function of A^qp for tt^’s in 200 GeV Au+Au and Cu+Cu collisions 
measured by PHENIX and charged hadrons in 2.76 TeV Pb+Pb collisions measured by ALICE. 


transverse to the beam axis satisfies the condition of d < The number of quark 

participants as a function of the number of nucleon participants is nonlinear, especially for low 
values of Npart- In Fig. El the similar feature as the previous plot is seen. Since the Apart is 
a factor of 2-3 higher than Aqp, all the points are systematically moved to the right. Finally, 
we plotted the dp^/pT against the charged multiplicity, dN^h/dr], as shown in Fig. [71 In this 
plot, we added one 62 GeV Au+Au point which is 0-10% centrality. We expect that dNch/dr] 
well represents the energy density of the system. It is interesting to note that the most central 
Au+Au 200 GeV points tend to merge into the most central points of Pb+Pb collisions, while 
they deviate each other as going to lower dNch/dr], This systematic trend has not been found by 
looking at Raa’s which look similar across the systems. In order to cross-check this new result, 
we have performed a power-law fit to dp^/pT vs dNch/dr] points from 200 GeV Au+Au collisions, 
and compared the power with the result obtained from a different method [9]. We fitted the 
points of this work with dp^/p^ = l5{dNch _/assuming dNch/dr] oc(Apart)^'^® [IH], and 
obtained a as 0.64+0.07. Assuming the spectra shape follows a power-law with the power n, 
one can write the relation between dp^/pT and Raa as: 

5loss = dpT/pT = /3Apart“, Raa = (1 - = (1 - /3Apart“)’""" 

Following this relation, we obtained the power a. as 0.57+0.13 from the fit to the integrated Raa 
as a function of Apart in the literature [Sj. We therefore confirmed that the powers obtained by 
the two methods are consistent. 

4. Summary 

We presented PHENIX measurement of the fractional momentum loss ((5 pt/pt) of high px 
identified hadrons. By looking at the (Ipt/pt instead of Raa-, we found many interesting features. 
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Pb+Pb 2.76TeV, Ch. hadrons: PLB 696, 30 (2011), dN /dr|: PRL 106, 032301 (2011) 
Pb+Pb 2.76TeV, Ch. hadrons: PLB 720, 52 (2013), dhi''/dr|: PRL 106, 032301 (2011) 
Au+Au 200GeV,7t“: PRC 87, 034911(2013), dN /dri: pSc 71, 034908 (2005) 

Au+Au 62.4GeV,7t“: PRL 109,152301(2012), dft /dr|: Preliminary 

Cu+Cu 200GeV,;i“: PRL 101,162301(2008), dN^/dTi: PHOBOS, PRC 83, 024913 (2011) 
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Figure 7. dp^/pT as a function of dNch/dr] for vr^’s in Au+Au collisions at 200 GeV and 
62.4 GeV, and in Cu+Cu collisions at 200 GeV measured by PHENIX and charged hadrons in 
2.76 TeV Pb+Pb collisions measured by ALICE. 


The dp^r/pT of high px which are computed from 39 GeV Au+Au over to 2.76 TeV Pb+Pb 
are found to vary by a factor of six. We plotted the dp^/pT against several global variables, 
-Vpart, A^qp and dNch/dp. It was found that 200 GeV Au+Au points are merging into the central 
2.76 TeV Pb+Pb points when plotting dp^/pT against dNch/dr}. We performed a power-law fit 
to the dpT/pT vs dNch/dr], and obtained a power that is consistent with the one obtained from 
the ht to the integrated Raa- We are going to add points from other systems to systematically 
investigate the dp^/pT. 
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